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Abstract 

We study the radiative decay B —>■ jiui in the framework of QCD 
factorization. We demonstrate exphcitly that, in the heavy-quark 
hmit and at one-loop order in perturbation theory, the amphtude 
does factorize, i.e. that it can be written as a convolution of a 
perturbatively calculable hard-scattering amplitude with the (non- 
perturbative) light-cone distribution amplitude of the 5-meson. We 
evaluate the hard-scattering amplitude at one-loop order and ver- 
ify that the large logarithms are those expected from a study of the 
h ^ u transition in the Soft-Collinear Effective Theory. Assuming 
that this is also the case at higher orders, we resum the large loga- 
rithms and perform an exploratory phenomenological analysis. The 
questions addressed in this study are also relevant for the applica- 
tions of the QCD factorization formalism to two-body non-leptonic 
5-decays, in particular to the component of the amplitude arising 
from hard spectator interactions. 



1 Introduction 



The study of -B-decays is providing a wonderful opportunity to test and 
improve our understanding of the Standard Model of Particle Physics and 
of its limitations. In particular, the BaBar and Belle 5-factories [|l], 0, as 
well as other experiments, are providing us with an impressive amount of 
accurate experimental information about two-body non-leptonic 5-decays. 
Among the remarkable achievements is the increasingly precise determina- 
tion of sin (2/5) (where (3 is one of the angles of the unitarity triangle) from 
measurements of the mixing-induced CP-asymmetry in the golden mode, 
B J/ipKs 0. Unfortunately the precision with which we can determine 
the fundamental properties and parameters of the standard model (in par- 
ticular the CKM matrix elements) from the measured branching ratios and 
asymmetries of other non-leptonic channels is severely limited by our inabil- 
ity to control non-perturbative QCD effects. Further progress in overcoming 
this limitation is urgently needed if we are to be able to exploit effectively 
the wealth of experimental data for fundamental physics. 

An important step towards the control of non-perturbative QCD effects 
in two-body non-leptonic 5-decays has been the recent discovery that in the 
heavy-quark limit, oo (where m;, is the mass of the 6-quark), hard and 

soft physics can be separated (factorized) Within this factorization 

framework the amplitudes are expressed as convolutions of perturbatively 
calculable hard-scattering kernels and universal non-perturbative quantities 
(light-cone distribution amplitudes and semi-leptonic form factors). When 
the decay products are light mesons Mi_2, such as in S tttt or _B ^ 
T^K decays, the factorization formulae take the generic form, represented 
schematically in fig.|l|, 

= y^Ff^^'\ml) [ duTl^{u)^MAu) + (Mi^Ma) 
Jo 

+ / didudvTl\i,u,v)^B{O^MAv)^MM- (1) 
Jo 

Here Oi is one of the operators of the effective weak Hamiltonian, p^^'^^^''^ 
denotes the B — >• Mi 2 form factors, and $x(w) is the light-cone distribu- 
tion amplitude for the quark-antiquark Fock state of meson X. Tl-{u) and 
T/^(^,M,t>) are hard-scattering functions, which are perturbatively calcula- 
ble. The hard-scattering kernels and light-cone distribution amplitudes de- 
pend on a factorization scale and scheme; this dependence is suppressed in 
the notation of eq. (|l]). Finally, mi 2 denote the light meson masses. In phe- 
nomenological applications of the factorization framework (see, for example 
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Figure 1: Graphical representation of the factorization formula. Only one of the 
two form-factor terms in eq. (|I]) is shown for simplicity. 

refs. 0, H), the hard-scattering kernels have been computed to 0{as). For 
the Tj^'s, for which the leading contribution is generally of order a^, this 
requires the evaluation of one- loop diagrams, whereas for the T/^'s, which 
start at 0{as), one only needs to evaluate tree- level diagrams. An important 
difference in the two cases is that the hard scale in the T^^j's is rrih, whereas 
in the T/^'s it is y^mftAqcD- The hard spectator interaction term in eq. (0), 
i.e. the component containing the T/^'s, therefore depends on three scales, 
''T^b, V^t^Aqcd and Aqcd, and it is important to understand the structure 
of the higher order corrections. Note that since we work at leading twist we 
do not distinguish between the quark mass rrib, the meson mass Mb or any 
scale which differs from these by O(Aqcd)- 

In this paper we study factorization, and in particular the one-loop correc- 
tions, in a simpler process, the radiative decay B — > 'jii'e, where the photon 
is hard. This process is interesting in itself, and also has the key features 
of the T/^ term in eq. (|l]), in particular the dependence on the three scales. 
The only hadron in the decay is the 5-meson, making it easier to focus on 
the properties of $b and on issues related to the presence of the three scales. 
The questions which we investigate include: 

1. Does factorization hold at one-loop order? 

This has not yet been directly verified for the T/^ term in eq. (|ID. In- 
deed it has been suggested that it is necessary to replace the light-cone 
distribution amplitude of the 5-meson by a wave function which also 
depends on the transverse components of momentum [0. 

With the definition of the S-meson's light-cone distribution amplitude 
given in eq. (H) below 0, Q , we verify that factorization does hold at 
one-loop order and that there is no need to introduce a dependence on 
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transverse components of momentum. We evaluate the hard-scattering 
kernel explicitly. 



2. How can one obtain information about the light-cone distribution am- 
plitude, ^B, non-perturbatively, for example from lattice simulations? 

For light mesons such as the pion, moments of the distribution ampli- 
tude are given in terms of matrix elements of local operators, and can 
be determined using non-perturbative methods such as lattice simula- 
tions or QCD sum rules. For the 5-meson the situation is different. 
The presence of ultraviolet divergences in the integral over a light-cone 
component of momentum implies that the corresponding matrix ele- 
ments of local operators do not give information useful for the decay 
amplitude. The distribution amplitude will therefore have to be deter- 
mined by evaluating the matrix elements of non-local operators, which 
is possible in principle, but considerably more difficult. Alternatively, 
one may attempt to determine the distribution amplitude from exper- 
imental measurements of processes such as the one being studied in 
this paper. The additional ultraviolet divergences also complicate the 
dependence of the distribution amplitude on the factorization scale (i.e. 
the evolution). The presence of such additional ultraviolet divergences 
has been stressed previously in ref. 0, and questions concerning the 
evolution of the distribution amplitude were discussed in ref. [0. 

3. How can one resum the large logarithms which appear in one- loop 
perturbation theory? 

In particular one has the Sudakov double logarithms associated with 
the heavy-to-light {b —* u) decay. Large logarithms appear both in 
the distribution amplitude and in the hard-scattering kernel. In this 
context we find the formulation of the soft-collinear effective theory 
(SCET) |]TU|, |TT| very helpful, and we present the resummed expressions 



in sec.|^. We stress however, that, whereas we have performed the one- 
loop calculation explicitly, the validity of the all-orders resummation in 
the present context still requires a formal demonstration. 

For some other processes, for example for the semi-leptonic B ^ n form 
factor at large momentum transfer, the lowest-order contribution in 
is singular at low momenta. In the Factorization approach of refs. [§, |^ 
the amplitudes for these processes are considered to be uncalculable. In 
the pQCD approach ||12| Sudakov effects are invoked to regulate these 



singularities and we have criticized the reliability of this procedure in 
ref. [0 (for a study of similar problems in higher-twist contributions to 



3 



the pion's electromagnetic form factor, see ref. ||14[)- For the B — > 'yiui 
decay the situation is different, the lowest order contribution has no 
singularity and there is no enhancement from non-perturbative regions 
of phase-space. 

4. What is a suitable choice of factorization scale fip"^ 

We argue in the following that it is (of order) ^m^AQCD- 

5. Should the light-cone distribution amplitude be defined in QCD or in 
the heavy-quark effective theory (HQET) ? 

The choice is a pure matter of convenience, since the decay amplitudes 
are independent of it. A different definition of the light-cone distribu- 
tion amplitude is compensated by a different hard-scattering kernel, in 
order for their convolution to remain identical in all cases. To connect 
our analysis with the results from the SCET - where heavy quarks are 
treated as in HQET - we find it convenient to define the distribution 
amplitude in the heavy-quark effective theory. In addition, since the 
factorization scale is much lower than rrib, it is natural to treat the 
light-cone wave function in the framework of the HQET. 

6. Are the one- loop results large? Does the resummation mentioned in 
item ^ make a large difference? 

We find that the one-loop corrections are typically of the order of a few 
times 10% and are sensitive to the choice of distribution amplitude. 
In particular, the corrections depend sensitively on the two parameters 
and A^'* defined in eq. (^8]). Resummation of the large logarithms 
typically changes the results by up to 30% or so of the one-loop correc- 
tion. 

In this study, we neglect terms which are suppressed by a power of m?,. In a 
number of the important two-body non-leptonic decay channels the higher- 
twist terms are "chirally enhanced" or have larger CKM- matrix elements and 
so can give significant contributions. The mass singularities in general do 
not factorize for the higher-twist terms, so that the Factorization framework 
described above has to be extended (for recent progress see ref. [|15[)- Here 
we focus on the higher-order perturbative corrections at leading twist. 

Some of the above questions have been investigated in previous studies, 
from which we have benefited considerably. Korchemsky, Pirjol and Yan 
performed a detailed study of this decay process in ref. [0 and concluded 
that a consistent factorization formula can only be derived at one-loop order 
if the hard- scattering kernel and the 5-meson distribution amplitude include 
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a dependence on the transverse momenta. We disagree with this conclusion. 
In sec. ^ we perform the matching exphcitly, with the distribution amphtude 
defined as in refs. @, §|, ^ Q from the non-local matrix element: 

= I dz.e^'-^- {0\M^)[z,OMO)\B)l^^,^^, , (2) 

where [z, 0] denotes a path-ordered exponential. (Here and for the remainder 
of the paper we find it convenient to write the label B on the distribution am- 
plitude as a superscript rather than a subscript.) We find that it is consistent 
for the wave function ($^) to depend on a single component of momentum 
/c+, and that the amplitude of the decay process can be expressed as a con- 
volution over the single variable 

^<p{K)T,.{h) . (3) 

Of course does depend on the details of the internal dynamics of the 
S-meson in a complicated way, which includes a dependence on the trans- 
verse momentum of its constituents. This does not matter however ($-^ is a 
quantity which must be determined non-perturbatively in any case). The key 
point is that the hard-scattering kernel only depends on We stress that, 
in general, fc+ should not be identified with a component of the momentum 
of any particular constituent of the 5-meson is defined through eq. (0) ). 

The remainder of the paper is organized as follows. In the next section 
we briefiy discuss the kinematics for the decay B — > We perform the 

matching, establish factorization and calculate the hard-scattering ampli- 
tudes at tree level and at one-loop order in secs.^ and || respectively. In 
sec. ^ we use the SCET to resum the large logarithms in the hard-scattering 
amplitude. We perform a brief phenomenological study in sec. | and in sec. ^ 
we present our conclusions and discuss open questions. 

2 Kinematics of the B — > ^^i^t Decay 

The B — > ^Ivi decay is illustrated schematically in fig. |^. The momentum 
of the i?-meson is denoted by = Mbv'^, where the four-velocity v satisfies 
= 1. In the following, unless otherwise stated, we will work in the rest- 
frame of the S-meson, so that v = (1,0). The momentum and polarization 

""^At leading twist, two distribution amplitudes can actually be defined from eq. (|^) for 
the _B-meson (see eq. (^) ). However, only one contributes to the B jii^i decay, and 
we shall call it the _B-meson distribution amplitude (or light-cone wave function) in the 
remainder of this paper. 



5 



B 




Figure 2: Schematic representation of the B '^Cmi decay. The momentum 
of the 5-meson is denoted by p and that of the photon (7) by q. 



vector of the photon are denoted by q and e* respectively. The energy of the 
photon is given by 

We consider decays in which the energy of both the photon and the lepton 
pair is large, of order M^, and we neglect mi /Mb, where mi is the mass of 
the lepton. 

The hadronic matrix element for the decay B — > jiu^ can be written in 
terms of two form factors Fy and Fa'- 

^ {^{e*,q)\u-f,{l-^,)b\B{p))= (5) 



e^,p,e*''v''q''Fv{E^) + i[e*^{v ■ q) - q^{v ■ . 

For light-cone dominated processes, such as the B — > ■jiui decay being 
studied in this paper, it is convenient to introduce the light-cone coordinates 
I — l±), defined by 

i± = ^^, l±^{h,k), l'^2l+l_-Q. (6) 

In the rest-frame of the i?-meson, we choose the photon's momentum to be 
in the minus direction so that: 

p = {Mb/V2, Mb/V2, 0^) and q = (0, g_, 0^) . (7) 
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We define the light-cone distribution amphtude of a state H which con- 
tains the 6-quark by 



dz. 



{0\Mz)[z,0%mH)\ 



z+,z±=0 ' 



where u,b are the quark fields and a,P are spinor labels. [z,0] denotes the 
path-ordered exponential V exp[—igs dz^A^{z)] and gs is the strong cou- 
pling constant. In the following we will use the terms light- cone distribution 
amplitude and light-cone wave function interchangeably. In both cases the 
terms denote <I'^^(A;+) as defined in eq. (|^). 

u current, 



is defined to be the matrix element of the weak b 



{^ie*,q)\u^,Lb\H) 



(9) 



where 'j^l = 7/^(1 — 75)- The question which we investigate in this paper 
is whether, up to one-loop order in perturbation theory, the matrix element 
can be written in the factorized form 



dkj^ 



(10) 



where the hard-scattering amplitude T does not depend on the external state 
[H) and is a function of hard scales only. We also evaluate T up to one-loop 
order in perturbation theory. 

Of course we are actually interested in the case where H is the 5-meson. 
However, in the perturbative evaluation of T we exploit the independence of 
T on the external state and choose a convenient partonic state H (in most 
cases we will take H to consist of a 6-quark and a w-antiquark). It is for this 
reason that we introduce the notation with the general superscript H. 

The matrix element on the right-hand side of eq. (|^) can be written in 
terms of two light-cone distribution amplitudes and and following 
ref. we write Q 



M{z) = {0\upiz)[z,0%{0)\B{p)) 



—I- 



JbM^ 



B 



1 + 



2$ 



B 



t 



75 



(11: 



a(3 



V ■ z and and depend on the coordinates z (they can be 

H= 



where t 

obtained from ^^^^ by Fourier transforming and projecting over the Dirac 



^Since in most of this paper we work in momentum space, we introduce a tilde on M, T 
and in coordinate space, and remove the tilde in momentum space. 
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(a) (b) 



Figure 3: Tree- level diagrams for the B 'yii'i decay with the two-parton 
external state \bu). 



structure). We will convolute this matrix element with a hard-scattering 
amplitude T{z): 

I d'zM{z)f{z) = I 0^J{k) J d'z/'-M{z). (12) 

If T{k) depends on a single component {k^ in our case), we can integrate 
the remaining components, which corresponds to setting z^ and Zx_ to in 
the matrix element M (the corresponding process is therefore a light-cone 
dominated one). We will see explicitly in the following sections that, at 
leading twist, it is only the function $:^(A;_|_) which contributes to the form 
factors Fv and Fa, and that Fy = Fa- 

It may be helpful to stress the distinction between A;_|_ and the kinematical 
variables of the initial state. The light-cone wave function depends on 
the latter in a complicated and non-perturbative way. An important goal 
of our investigation is whether it is sufficient to introduce a light-cone wave 
function as in eq. (||) depending only on the + component of k (as well as on 
all components of the momenta in H) or whether a generalization of eq. (P) is 
necessary. We will see that, up to one- loop order at least, eq. (H) is sufficient. 

3 Factorization at Tree Level 

In this section we study factorization for the process B ^ivi at tree level. 
This is straightforward and the result is well-known []^, |TB[. However, the 
calculation is instructive and so we exhibit the ingredients systematically. 
This will help make our presentation of the one-loop calculation in sec. ^ 



8 



below clearer. Of course many features of the generic calculation are absent 
at tree level. We do not encounter mass singularities so that factorization 
at tree level is guaranteed, and we do not have ultraviolet divergences and 
hence do not have to choose a factorization scale. These questions will arise 
at one-loop order in sec. |[ 

We wish to evaluate the hard-scattering kernel at tree level. For the 
factorization formalism to be applicable, the hard-scattering kernel must be 
independent of the infra-red effects contained in the initial state and we are 
therefore free to choose this state in any convenient manner. In this section 
we take the initial state to be a light antiquark (u) with momentum k (with 
the components of k of order Aqcd) and a 6-quark of momentum p — k. The 
tree-level result for the hard-scattering amplitude is given in eq. (|T^ below. 
It is also instructive to see how the same hard-scattering amplitude arises 
with a different initial state. In sec. |3.1| we illustrate this by obtaining the 
result of eq. (^7^ with a three-body, bug (where g represents a gluon), initial 
state. 

The two diagrams which contribute to the form factors at tree-level are 
represented in Fig. ^. In diagram |^(a) the (internal) 6-quark propagator is of 
0(l/mb) whereas in diagram |^(b) the internal propagator, which is now the 
propagator of the w-quark, is of 0(1/Aqcd)- Thus diagram (a) is suppressed 
by one power of the heavy-quark mass relative to diagram (b) and will be 
neglected in the following. Also, in the evaluation of one-loop graphs in 
sec. ^ below, we will only need to consider the diagrams in which the photon 
is radiated from a light-quark. 

In order to evaluate the hard-scattering amplitude we proceed in the 
standard way: 

1. compute the matrix element corresponding to the diagram of fig ^(b); 

2. evaluate the light-cone wave function at tree level; 

3. combine the two to deduce the hard-scattering amplitude. 
We now carry out each of these three steps in turn. 

The Matrix Element at Tree-Level: At leading order in l/rrif,, dia- 
gram ^(b) gives the following contributions to the matrix element 

Fjfi)"^ ^ {^{e*,q)\u^,Lb\b'ip-k)u%k)) 

= {v^{k)f^^,,u'{p-k)}, (13) 

where u and v are the spinor wave functions of the b and u quarks respectively 
and s and S are spin labels. e„ is the electric charge of a w-quark. 
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The Light-Cone Wave Function at Tree Level: We define the hght- 
cone wave function, of the initial state consisting of a 6-quark and 

a M-antiquark as in eq. (|^): 

^%{k^) = I dz. e^^+^- {0\M^)[z,0]b^{0)\b'{p~k)u\k))\^^^^^^^ . (14) 

depends on the initial state (for example in this case it depends on k), 
but we leave this dependence implicit. At tree level the wave function is 
readily found to be 

$l°j'"(fc+) = 27r5(A;+ - ~k^) v^^k) ui{p - k) , (15) 

where the superscript (0) denotes tree level. 



Matching and the Determination of the Hard-Scattering Ampli- 
tude: Writing the matrix element in the factorized form: 

i^r"= j^-^^T^h)Tn{h). (16) 

and taking the expressions eqs. ([13D and ([15|) 

respectively, we obtain the hard-scattering amplitude at tree level, 



Tfi{K) = --\[ni,LW. (17) 



As expected the hard-scattering amplitude depends only on k^ (through the 
hard variable q ■ k). 

3.1 Determining Tp'^ik^) with a 3-Body Initial State 

For the factorization formalism to be valid the hard-scattering amplitude has 
to be independent of the external state. For example it should be independent 
of the choice of external momentum k (which we verify at one-loop order in 
sec. ^). In this section we verify the independence of T*^"^ from the external 
state in an instructive example. We take a three-body initial state (see fig. §) 
consisting of a w-quark with momentum fc, a gluon with momentum / and a 
6-quark with momentum p — k — I. Now the matrix element at tree- level is 

pmm, ^ {',ie\g)\uj,r.b\b{p-k-l)u(k)g{E"'>,l)) 
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Figure 4: Tree-level diagram for the B 'yii'e decay with the three-parton 
external state \bug). Only the leading twist contribution is shown. 



where gs is the strong coupling constant and e^^^ is the polarisation vector of 
the gluon. 

The lowest-order contribution to the light-cone wave function of the three- 
body initial state is of 0{gs) and corresponds to the external gluon being 
annihilated by a gluon field present in the path-ordered exponential [^,0]. 
The tree-level term in the wave function is 

= -^gsv',{k)u'^{p-k-l)e^f{l) 

X J dz^j c/a2_e-^('=++"'+-^+)^- (19) 

= (27r) [5{k+ + /+ - h) - 6{k^ - h)] ef{l) v^i^k) ul{p - k - I) . (20) 
'+ 

The integral over a in eq. (0) is over the position of the gluon field which 
annihilates the incoming gluon (this position is taken to be az). 
It is now straightforward to verify that 

^(0)...^ (21) 

where T^°-*(A;+) is given in eq. (0). Thus we have verified that, as required, 
the same hard scattering amplitude is obtained also with the three-particle 
initial state. 

3.2 B -fiui at Tree Level 

Now that we have calculated the hard-scattering amplitude at tree level, we 
can express the form-factors defined in eq. (^) in terms of the light-cone wave 
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function of the i?- meson. In this way we obtain the standard result: 

= Fv = r + 0{a., 1/Ms), (22) 

Zy/Zh^ Jo 

where Qu = —2/3 is the charge of the u antiquark in units of the proton's 
charge. At this order the two form-factors are equal and only depend on 
the first inverse moment of This is similar to the corresponding 

calculations of amplitudes of two-body non-leptonic i?-decays. The leading 
contribution to the term on the second line of eq. (P also depends on the 
first inverse moment of $5. 



4 Factorization at One-Loop Order 

In this section we establish factorization for B '^ift decays at one loop 
in perturbation theory and calculate the hard-scattering amplitude at this 
order. We expand the matrix element, wave function and hard-scattering 
amplitude in perturbation theory, so that the factorization formula takes the 
schematic form, 

= F^o)^ + F«^ + --- = <l>^®T (23) 
= [$(0) ^ ® T(°)] + [$(°) ^ ® T(i) + $(1) ^ ® T(°)] + . . . 

where ® denotes the convolution, and the superscripts indicate the power of 
as- The hard-scattering kernels T^"^ contain only hard scales, whereas the 
distribution amplitudes absorb all the soft effects. 

At tree level we did not have to specify precisely the criterion used to sep- 
arate "soft" and "hard" scales. As discussed in the introduction, B — >• 7^//^ 
is a three scale process, with a large scale (mf,), a small scale (Aqcd) and 

an intermediate scale {^2q ■ k where k = (A;+, 0, 0_l)). We assume that large 

values of /c+ are damped by the hadronic dynamics so that fc+ is predomi- 
nantly of O(Aqcd) and the convolution is dominated by the region in which 
q ■ k = 0{mbAQCB), and that mftAQco is a sufficiently large scale for pertur- 
bation theory to be applicable. We will argue below that it is convenient to 
choose the factorization scale /x/ to be of 0{^ym^,AQCB)■ In this section we 
evaluate T^^^ and demonstrate that it is free of any infrared effects. T^^^ does 
however, contain large logarithms which will have to be resummed. This will 
be discussed in sec. |. 

As in sec. ^, we evaluate T^^^ by taking H = \b^(p — k)u^(k)). The 
components of k are taken to be of O(Aqcd), but the result for T^^^ does not 
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depend on the precise choice. T^^^ is obtained by using eq. p3| ) 

$(0)^®T« = F^i)^-${i)^®r(°), (24) 

so that, at one-loop order we need to evaluate both <|>(^)-'^ and Fjp^ . There 
are mass singularities present in both the terms on the right-hand side of 



eq. (|2^) , but we shall demonstrate that they cancel in the difference. Indeed, 
as will become clearer below, it is convenient to consider the difference in 
eq. (^) diagram by diagram Q The mass singularities cancel diagram by di- 
agram, and it is of course necessary to regulate them in the same way in both 
F^^^ and ®T'^^\ Unless explicitly stated to the contrary, we regulate 

the coUinear divergences by giving the u antiquark a mass (fc^ = m^) and 
the infrared divergences by giving the gluon a mass A. The hard-scattering 
amplitude is independent of m and A. The ultraviolet divergences are regu- 
lated by dimensional regularization (we work in d = A — e dimensions) and 
we use the MS renormalization scheme, by redefining /i^ — >■ fi^e'^'^ /Air, where 
je is Euler-Mascheroni constant, and subtracting divergences proportional 
to powers of N^^ = 2/e. (At one loop, if only single poles are present, this 
involves the subtraction of terms of the form 2/e — je + log(4vr).) 

Eq. (|2^), together with the tree-level expressions for the light-cone wave 
function for the bu initial state in eq. (0) and the hard-scattering amplitude 
in eq. (^^, gives the following expression for T^^^ 

v^{k)T^^^u\p -k)= F^^^'- - I '^^^T^k^)v^\ril,LW . (25) 
We now evaluate the contributions to each of the two terms on the right- 



hand side of eq. (p5D from each of the diagrams. Since all the calculations are 
performed with an external state of a 6 quark and a u antiquark, for com- 
pactness of notation we will omit the hu label on F and $ in the remainder 
of this section. 

4.1 Electromagnetic Vertex 

The contribution to F^^'* from the vertex correction to the electromagnetic 
vertex (the diagram in fig. ^(a) ) is 



cm 



X 



^ 2q_k+ J {2tiY 

v^{k)Yilt - ^ -^)ipii,Lu'{p - k) 

{P + ie){{k + If - m2 + ie){{q - k - If - m? + ie) 



(26) 



^We perform the calculations in the Feynman Gauge. 
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(a) (b) 



Figure 5: Electromagnetic vertex: contributions to the form factors (left) and 
^(i)em ^ rp{o) (^Yight). The dashed line in the right-hand diagram represents 
the path-ordered exponential. 



where we have neglected terms which are suppressed by AqcD/^b- Evaluat- 
ing the integral, we obtain 



p(l)em 



a.Cf 



Att 2q_k 



- log — ^ + 2 log ^ 



(27) 

where fiR is the renormalization scale (since the electromagnetic current does 
not get renormalized, the dependence on fiji cancels with the corresponding 
contributions from the wave function renormalization) |^. If we take the small 
mass limit by decreasing k_ and k± while keeping k^ fixed (and of O(Aqcd) )? 
then the mass singularities of the form log(m^) in eq. (|27| ) come from the 
coUinear region of phase-space, /+ = 0(A°), /_ = O(A^) and l± = 0(A) with 
A ^0. 

The corresponding contribution to the distribution amplitude is: 



$(l)em 



daz^ 



{2i 



{2nY 



-iz- {k+l—a /)_!_ 



7+ 



-iCpg. 



dH 1 



(27r)^/+ 
x[r7+(/ + ^)]/3nf 



{2Ti)[5{k+-k 
1 



- 5{l 



sl_ 



+ k+- ~k+)] 



[{k + - m^][P 



(29) 



■*In table 1. of ref.Q] there is an additional —1 in the square parentheses. We have 
traced this discrepancy to the expression for in eq. (A9) of ref. M. 
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a is the position of the gluon field in the path-ordered exponential (drawn 
as a dashed line in Fig. ^(b)) used to construct the gluon propagator. After 
the integration over a, followed by that over z_, we obtain eq. This 
equation underlines again that k+ cannot be identified with and is not a 
priori of O(Aqcd) (it is fc+ which is taken to be of O(Aqcd))- However, the 
convolution with T^^^ (which is proportional to l/A;^) damps large values of 
A;+ and is indeed dominated by the region in which is of O(Aqcd)- This is 
also true for the other one-loop diagrams, as can be seen from the expressions 
given in the following sections. 

Using the expression for <|)(^)'^™ in eq. (pS]), the contribution to ®T'^^^ 
is found to be: 



2q^k+ J {2tiY 
v%k)Y{t^ +f)filpil^.LU^{.P - k) 
{P + ie){{k + lf-m? + ie){-2q^{k+ + l+))' ^ ' 

Comparing eqs. (^) and (^) we find that, as expected, they are related 
by the replacement of the internal light-quark propagator by the eikonal 
approximation: 

{q~k~iy-m^ -2q_{k+ + l+)' ^ ' 

The diagram in fig.^(b) represents <|)(^)''™ (g) T^'^\ with the dashed lines rep- 
resenting the eikonal propagators. 

In the coUinear region of phase space, the replacement in eq. (|3T1) could 
be made also in f'^''^™ and therefore we would expect that the collinear 
singularities are the same in pj^^'^'^ and T'^^\ and this is indeed the 

case. Evaluating the integral in eq. ( pUD we obtain the result 

asCp 



^(l)em ^ rjn{0) 



2N^y - 2 log 



m 



+ 4 



(32) 



4:71 2q_k 

The integral has ultraviolet divergences which again we regulate using di- 
mensional regularization and renormalize using the MS scheme with a renor 
malization scale ^p. iip is the natural choice for the factorization scale. 

We see that, as anticipated, the mass singularities in F^^^*^™ and 
T*^°^ are the same, and from their difference we obtain the following contri- 
bution to the hard-scattering kernel: 



(33) 




15 



Figure 6: Wave-function renormalization diagram for the internal light-quark 
line. 



4.2 Wave Function Renormalization 

The contribution of the diagram in Fig. to F^^^ can be readily evaluated: 



p(l) wfc 



4:71 2q_k^ 



-Nr + log 



2k ■ q 



1 



(34) 



There is no corresponding contribution to The gluon propagator in 
this case is (0|yl+(az_)yl+(/5z_)|0), where az_ and I3z- are two points on the 
line between the origin and z (a, (5 are integrated between and 1, respecting 
the path-ordering). In the Feynman gauge this propagator vanishes. 

The contribution of Fig. |^ to the hard-scattering kernel is therefore: 



(35) 



which is free of mass singularities. 

We now turn to the wave function renormalization of the external quark 
fields. We start by considering the wave function renormalization of the u 
antiquark, 




^(i){iwf ^ 2.(0) 



i[Zf-l](/.«)jF(°) 



(36) 
(37) 



The wave function renormalization constant Z2 for a quark with momentum 
p and mass m in QCD is defined here in terms of the one-particle irreducible 
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graphs S by 



l + i 



(38) 



Since infrared effects cancel in the difference of the terms in eqs. (^) and 
(|37|), it is convenient in this case to regulate the infrared divergences by 
choosing the u antiquark to be off-shell, k"^ ^ m^. Then 



2 ' 





m 




+ log--2 


\ ) 





(39) 



The corresponding contribution to the hard-scattering kernel is 




(40) 



For the renormalization of the external 6-quark field, the contribution to 
^ii)Qw{ ^ rp{o) corresponds to the HQET. The superscript Q denotes heavy- 
quark, as defined in the HQET, and is used to distinguish it from the 6-quark 
in QCD. We have 



^(i)Qwf^y(o) ^ <!i[Z2«-i](^^)l$(o)®T(°) 



(41) 
(42) 



where Z2 is given by eq. (|39D which we now rewrite in the form (we take 
= ml) 



i|z5-i]W) = ^ 



2 log: 



+ 31og— -2 



(43) 



For a quark Q with residual momentum k in the HQET we define the wave 
function renormalization constant through 



l+i 



d{v ■ k) 



v-k=Q 



and find 



iVf^ -21oj 



2v k 



(44) 



(45) 
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After renormalization, the contribution to the hard-scattering kernel is 




(46) 



We end this subsection by using dimensional regularization to regulate 
the infrared divergences, instead of taking the external quarks to be off-shell. 
This provides further verification of the independence of the hard-scattering 
amplitude from the regularization, and allows us to compare our results with 
earlier studies. Using dimensional regularization in QCD at the scale fi we 
find: 



An 



-In, 

2 ' 



uv 



3 m 

+ 77 log — 
2 fi' 



IR 



(47) 



for a quark q of mass mq {q = uoi h) with N^^' = 2/5 (recall that a rescaling of 
/i has been implicitly performed and cancelled the contribution proportional 
to —7b + log(47r)). The result in eq. (|i7|) agrees with eq. (29) of ref. |]1T|, but 
disagrees with that in table 1 of ref. . For the heavy quark effective theory 
we find 

\[Z^-lW = '^[Nr-Nl^] . (48) 

To compute the contribution of the 6-quark wave function renormalization 
graph to Th in dimensional regularization, we take the difference of eqs. (|47|) 
(with the quark mass set to mi,) and (^8]). Comparing the result with 
eq. (^61), which is the same contribution evaluated in the off-shell regulariza- 
tion scheme, we see that it is indeed the same for /x = /i^? = /i/j Repeating 
the analysis for the wave function renormalization graph for the light quark, 
we find vanishing contributions in both cases. Thus, as expected, we obtain 
the same result for the hard-scattering kernel for both methods of regulating 
the mass singularities. 
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Figure 7: Diagram representing the one-loop correction to the weak vertex. 



4.3 Weak vertex 

We now turn to the correction to the weak vertex illustrated in the diagram 
in fig. 1^. Evaluating this diagram we find 



p(l) wk 



47r ^ 



X y . y . 2y 

log - - log - 

X z 



/i^ x-y 



(49) 



+2 log ^ log ^ + 2 log ^ + log' - + 2Li2 f 1 - ^) + 4Li2 

y z z y \ xJ \ yj 



+ . 



with 



X 



m 



b ) 



y 



and 



2(g ■ k) 



(50) 



The ellipses denote terms which are not relevant for the discussion, since 
their Dirac structure F is such that they vanish when convoluted with the 
distribution amplitudes: Tr[(l + |^)75F] = Tr[(l + 1^)7+ 75F] = 0. 

In ref. [0 this diagram was evaluated in the HQET. The infrared be- 
haviour of eq. (ff3): 



log^ fc+ + 2 log k+ \og{V2E^) - 2 log k^ 



(51) 



agrees with that in ref. 

^In the case of dimensional regularization, there is only one scale fj, for both ultraviolet 
and infrared divergences. We need therefore to set the factorization and the renormal- 
ization scales to the same value in order to compare the hard-scattering kernel with the 
result obtained in another scheme. 
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The corresponding contribution to ® T^^^ is 

agCp 



^(l)wk^y(O) 



X 



An 



2iVf^log 



2q_k 



(52) 



2 log 



2 



+ . . . 



We recall at this point that we are using the modified minimal subtraction 
scheme (MS). We have therefore redefined /z^ /i^e'^^/47r, and will subtract 
the divergences proportional to powers of N^^ = 2/e. 

The contribution to the hard-scattering kernel contains double logarithms 
and is given by: 




(53) 



4.4 Box diagram 

The final one-loop graph is the box diagram in fig. |^. It has been argued in 
ref. that the structure of this graph requires the introduction of a trans- 
verse momentum dependence into the 5-meson wave function, and hence 
invalidates the factorization formalism encapsulated in eq. ([lOD . We disagree 
with this conclusion. We will see that the leading-twist contribution from the 
box diagram is a soft effect which is absorbed completely into the light-cone 
wave function (i.e. the contributions to F^p and ® T^^^ are the same). 
The question is not whether the box diagram exhibits a dependence on the 
transverse components of the external quarks (which in general it does), but 
whether this dependence matches the one induced by the distribution am- 
plitude. This is indeed the case and therefore we find no breakdown of the 
QCD factorization framework at one-loop for B —>■ jiue decays. 
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Figure 8: Box diagram. 



The contribution to Pji^^ from the box diagram in fig. ^ is 

[(/ + ky — m? + ie] [P + ie] [{q — k — ly — m? + ie] [{p — k — ly — ml + ie] 

We are interested in the leading twist contribution to Fj^\ which is of 
0{l/k^). At one- loop level the box diagram is the only one with no propaga- 
tor of order 1/Aqcd outside of the loop. The l/fc+ enhancement is therefore 
more difficult to obtain than for the other diagrams; it must come from sin- 
gular regions of phase space in the loop integral. We now investigate what 
these regions are. 

Consider the region of phase-space in which the components of / are of 
order ~ /;;„ ~ Aqcd- Power counting shows that this region gives a 
leading twist contribution. However in this region the 6-quark propagator 
can be approximated by that of the HQET and 

i-if^-l^ + m i ^^^^ 



{q-k-lf-m? + ie -2g_(A;+ + /+) ' 

i.e. the light-quark propagator joining the electromagnetic and weak ver- 
tices can be replaced by its eikonal approximation. We therefore find that 
the contribution from this soft region is equal to the corresponding term in 
$(i)box ^ 2^(0) _ jn(igg(i obtain at leading twist 

^a)box _ ^(l)box^y(O)^ (56) 
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so that the corresponding contribution to T^-^^ is zero, 



r«'^°^ = o. (57) 

We conclude this subsection by briefly illustrating the suppression of con- 
tributions from regions of phase-space which would lead to a non-zero contri- 
bution to T^^^^°^. As an example consider the contribution from the collinear 
region in which /_ ~ 0{mf,), /+ ~ O^AQ^j^/nib) and l± ~ O(Aqcd)- In this 
case the phase-space is of 0(Aqqj-,) and the four factors in the denominator 
in eq. (0) are of order AqcD^^fe, ^qcd; ^qcd?^6 and ml respectively. This re- 
gion therefore does not give a contribution of 0(1/Aqcd) (in addition there is 
a further suppression factor of Aqcd from the numerator). A detailed study 
of fIi^^^°^ and $(^)''°^ (g) T^"-' shows that they are equal at leading twist. 



4.5 The Hard-Scattering Kernel 



We obtain the complete one-loop hard-scattering kernel by summing the 
contributions from the electromagnetic vertex, eq. (|33D, the wave function 
renormalization of the collinear light quark, eq. (^), the wave function renor- 
malization of the external w-antiquark and the 6-quark, eqs. (^0|) and (|46|) 
respectively, the weak vertex, eq. (|53|) and the box diagram eq. ([571) . We 
now convolute the result with the wave function of the 5-meson, using the 
decomposition in eq. (pAf), obtaining the leading twist one-loop expression 
for the form factors: 



/ 



(5^ 



where 



T{k+,E^-^iF) 



JeMbQu 1 
2V2E^ k+ 



(y. C 

1 + -^^K{k+,E^]HF) 



(59) 



and 



K{k+,E^;^F) 



log' 



2fc+g_ 1 



mt 



mt 



-log^^ + -log^ + 21og^log 



-2 log 



+2Li2 



V2q- 



V2q^ mb 



lot 



4Li, 1 



rrib 



V2q- 
rrib 

rrib 



V2q- 



^^F 



V2q- 
(60) 
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The dependence on renormalization scale /x/j has canceled as it must do. 

The formulae in eqs. (|58D - (|60D above give the one-loop factorized ex- 
pression for the form factors and Fy. The expression for T in eqs. ( |55| ) 
and ( |60D differs from the corresponding result in refs. [0, [1^. Of course T 
depends on the definition of the light-cone wave function and we present 
our definition in sec. |^. With this definition there is no need to introduce a 
dependence on transverse momenta, at least at one-loop order. We disagree 
also with the formula proposed in eq. (12) of ref. []16| . 

The large double and single logarithms present in eq. log^(A/2fc+/mfe) 
and log(A/2A;+/mb), need to be resummed and this is the subject of the next 
section. 



5 Resummation of Large Logarithms 

In order for the expressions for the hard-scattering amplitude to be use- 
ful phenomenologically, the large logarithms need to be resummed. This is 
true in particular for the Sudakov effects associated with the 6 — > m weak 
vertex. An elegant method for performing this resummation exploits the 
effective field theory describing the interaction of (infinitely) heavy quarks 
and (massless) coUinear quarks with soft and collinear gluons |Ty, |TT|. Al- 
ternatively one could perform the resummation by using the "Wilson-Line" 
formalism ^ . We now adapt the discussion of ref. |TI| to the B •yiug 

decay process we are considering. We should stress however, that, whereas 
the one-loop calculations in the preceding section were explicit, here we are 
assuming that we can incorporate our one-loop results into the general frame- 



work of refs. 1 10, 11 1. This needs to be verified explicitly. 



5.1 The Matrix Element 

In this section we consider the resummation of the large logarithms in the 
amplitude F^. We start by separating the total expression for T(i) into 
a contribution from the electromagnetic and weak currents, including the 
corresponding contributions due to wave function renormalization (since the 
contribution from the box diagram is zero we do not discuss it anymore in 
this section): 

2^(1) _ 2^(1) Jem _j_ 2^(1) Jw ^Q^^ 
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where 



em _j_ ^'T'(l) wfc _j_ T^(l) 

2 



uwi 



4:71 



3, 2k -q 9 



(62) 



and 

'IpWJw — 

V2g- 

+ log' 



T-i(l)wk _|_ ^T-i(l)wfc 1 T-i(l) 6wf 

2 



47r 



2^(0) 



2fc • 



- log —j^ + log ^ + - log2 ^ 



log + log 



r ^ i^B —j^ log —i=^ y 2 log 

V2A;_|_ \'2q_mb v2A;+ 




(63) 

rrib 
V2k+ 



+ 4Li2 1 - 



V2q-. 



From eq. (|62D we readily see that if we choose the factorization scale fip to 
be of 0{q ■ k) = 0{mbAqcD), then there are no large logarithms in T*^^^'^'=™. 
On the other hand, this is not the case for T^^^'^^. 

We focus therefore on the contribution from the weak vertex, which in- 
volves a transition of the 6-quark with momentum p — A; to a w-quark with 
momentum g — /c: 



{u{q-k)\uj^{l-j,)b\b{p-k)) 



(64) 



From eq. (^3]) we see that even at one-loop order the weak transition vertex 
has large double and single logarithms, and it is these logarithms which 
we attempt to resum using the techniques of the SCET |TD|, |rT|. Thus we 
are considering the set of diagrams, where an arbitrary number of gluons is 
exchanged between the incoming heavy quark and the collinear light quark 
that links the weak and electromagnetic vertices. The effective theory is 
designed to describe processes with heavy quarks and (almost light-like) light 
quarks. Fluctuations below some scale /x^ (<^ ml) are described in terms of 
the effective theory and contributions from higher momenta are obtained by 
perturbatively matching onto QCD and absorbed into coefficient functions. 
In addition to heavy quarks, the effective theory contains two kinds of light- 
quark and gluon fields. Using light-cone coordinates (writing an arbitrary 
momentum p as p = {p+,P-,p±)), the effective light-quark and gluon fields 
are either soft / oc mb(A^, A^, A^) or collinear / oc mb(A^, 1, A), with A a small 
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expansion parameter. In our case, the collinear quark carries momentum 
q — k, and thus A = 0(a/Aqcd7"^)- We are interested in the leading- twist 
contribution to the matrix element, and therefore in the SCET at leading 
order in A as discussed in ref. (for an extension to higher orders in A, see 
ref. m). 



The weak current in eq. (|6^) can be matched onto the operators Oi of the 
effective theory: 

uY{l-l5)b = J2c^{^^)o^{^^)■ (65) 

i 

The Ci are Wilson coefficients which describe the physics above fi and depend 
only on the hard scales such as rrih and g_. The dependence of the Wilson 
coefficients, Ci, on fi cancels the one which is implicit in the definition of the 
operators. The large logarithms contained in the coefficient functions can be 
resummed using renormalization group equations. As discussed below, this 
means that the Wilson coefficients will receive contributions from all orders 
of perturbation theory. 

For the process which we are considering, the B — > 'yiui decay, only two 
SCET operators are relevant O3 and Oe, which correspond to the coefficient 
functions C3 and Cq in eq. (27) of ref. |Tl|. These two operators have the same 
Wilson Coefficient C^fip) = C^^jjip) = Cq^hp) (in fact, only the difference 
O3 — Og is relevant for the left-handed current in this process). Thus we need 
the renormalization group equation for a single coefficient function. 



In Ref. |]Tl| one-loop results were obtained for the weak current. When 
convoluted with the lowest order wave function and rewritten in our notation, 
these results are: 

_p(0) I pJw ^ _p(0) I _p(l)wk , 1 _p(l)wfc , _p(l)fewf 

fl ' ^ fJ, ' l-L ' fl ' fl 



(66) 
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where 

Lf(M = -(iVf^)^-2iVf^log-g- (67) 

^/2k+ 4 

L^^M = ^f'' + 21og-^ (68) 
Lf(M = 2(iVf^)2 + 2iVf^ + 2iVf^log^^ (69) 

+ ^-^2(f^ + ^^-2(f^ + ^-Y 
^ (70) 

Lg {Lc) comes from the exchange of a soft (coUinear) gluon, and the super- 
scripts "wk" , "6wf " and "wfc" denote respectively the one-loop contributions 
from the weak vertex, the wave function renormalization for the heavy quark 
and (half of) that of the coUinear light-quark. One might have also expected 
two additional terms on the right-hand side of eq. (§^), LJ^'^ and L'^^. How- 
ever LJ^^ and L^"^ vanish in the effective theory because of the decoupling 
of soft gluons and coUinear quarks and of coUinear gluons and heavy quarks 
respectively. 

At this point we should point out that the result for L^™^ differs from 
that in eq. (35) of ref. [|1T| by a constant term. This is due to our choice for 



the definition of Z2 when the infra-red divergences are regulated by taking 
the quark off-shell (we define Z2 through the derivative of the self-energy 



diagram with respect to the external momentum as in eqs. (|38|) and (|4^) ). 
The difference cancels in the matching, as long as one consistently uses the 
same definition in the effective theory and QCD, leading to the same result 
for T(i). 

The matching of SCET onto QCD at the scale /i = rrih yields the Wilson 
coefficient at this scale: 

C{mh) = 1 — (71) 



X 



21og'^ + 2Li2 1-- +^ log^ + — + 6 

X \ x/ X — y X 12 



where x and y have been defined in eq. (^0]). 

We now consider the renormalization group equation for C (fi) and study 
its scale dependence. From the one-loop results in eqs. (§3)"(E0)' obtains 
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the evolution equation for the Wilson coefficient function [|Tl|] 

7(/i)C(/i) , 



where the first two terms of the anomalous dimension 7 are given by 



7lo 



iNLO 



TT 



log 



Air 



fx 



V2q. 



(72) 

(73) 
(74) 



with B = Ca(67/18 - 7i^/6) - 5iV//9. The value of B is deduced in ref. fTT 
by comparison with earlier work on inclusive B Xg'y and B X^ii^ 
decays |18|, |19|. The coefficient B can (and should) be checked directly by a 
two-loop computation within the SCET framework. 

The homogeneous nature of the evolution equation leads to the exponen- 
tiation of (Sudakov) logarithms: 

C(/i) = C{mb) exp 



asirrib) 



where 



Mr) 
Mr) 



1 + log r 



1 — r + r log r 



log^ r 



(75) 

(76) 
(77) 



■^-21og^' 
2 



logr 



2CfB 



[r — 1 — log r] 



with 



(7J 



r = asijx) /asirrib) 

(3o = HCa/S - 2Nf/3 and (3i = 34C1/3 - 10C^A^//3 - 2CFNf. 

In fig. H the Wilson coefficient C (/i) is plotted as a function of fj, (at next- 
to-leading order as in eq. (|75D ). At this order, C is also a function of the 
photon energy Ej, and for illustration we chose two different values of E^. 
We note that for of order nib Aqcb, C differs from 1 typically by 10-20%. 

The above equations allow us to resum the large logarithms present in 
T^^\ Before doing this however, we verify that our results at one- loop order 
agree with those of ref. Expanding the exponential in eq. (ffSD we obtain 

Cifip) = C{mb)x (79) 
I -^log- 

rrib 



1 + 



An 



+ 4 log log ■ 



rrib 



rrib 



5 log 



rrib 



0{al) 
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Figure 9: Value of the Wilson coefficient C at next-to-leading order as a 
function of /i, for two different values of the photon energy: = Mb/ 4: 
(dotted line) and = Mb/S (dashed line). 



We now need to combine the one-loop term in eq. ( ffOD with the one-loop 
expression for the coefficient function at the scale rrib in eq. (|7Tp and the re- 



maining one-loop terms in eq. (^), using the explicit expressions in eqs. 
(|7D|). The dependence on the factorization scale disappears as expected, and 
the result is identical to our result for fI°^ + F^^ , see eqs (|T|), (|3|), (|T|) 
and (M). 



5.2 The Distribution Amplitude and Hard-Scattering 
Kernel 

In the previous subsection, we have checked that the SCET correctly repro- 
duces the one-loop contribution of the weak current to the matrix element. 
In this section we investigate whether the contribution to the distribution 
amplitude can be understood in the effective theory as well. We shall see 
that this is indeed the case, allowing us to give an expression for the hard- 
scattering kernel in which the large Sudakov logarithms are resummed. 
Recall that the contribution to the distribution amphtude corresponds 
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to Feynman diagrams where the propagator of the internal hght quark is 
replaced by its eikonal approximation. In the framework of SCET, this ap- 
proximation is automatically performed when the interaction of soft gluons 
and coUinear quarks is considered. We expect therefore that the contribution 
to the distribution amplitude in our framework is equal to the exchange of a 
soft gluon in the SCET framework. 

When we compare our results for the contributions to ® T^^^ in 
eqs. (H) and (^) with those for L^™^ and LJ^ in eqs. (^) and we 
find explicitly that: 

$(l)Qwf ^ j.(0) ^ ^2^(0) ^^wf (gQ) 

and 

${l)wk ^ y(0) ^ ^2^{0)(^^) ^wk ^ ^g^) 

An 

for any value of the factorization scale. Both ^i^^)™*^^ and LJ^"^ are zero. The 
soft terms, the L^'s, do therefore correspond to the terms absorbed in the 
distribution amplitude and the hard-scattering kernel contains the remaining 
terms, i.e. the collinear contributions and the resummed logarithms in 
the Wilson coefficient C. 

At this point it may be instructive to consider the various scales which are 
present in the calculation, and we illustrate the corresponding contributions 
to the matrix element in Fig. Fig. ^(a) illustrates the situation for the 
electromagnetic vertex. The contributions from scales between Aqcd and 
y^mbAqcD are "soft" , in the sense that all components of momenta are small 
and the diagrams satisfy the SCET Feynman rules for soft gluons (although, 
as seen in sec. [4.1| , in the evaluation of diagrams the components are not 
always uniformly small). The mass singularities and remaining contributions 
from scales between Aqcd and \/ rrib Aqcb, are absorbed into the distribu- 
tion amplitude, so that hard-scattering kernel has no large logarithms for 
fip = O { nib Aqcb)- The situation is more complicated for the weak ver- 
tex in Fig. [Ty(b). The matrix element gets large logarithms from soft and 
collinear regions of phase space. The soft contributions come from scales 
between nib and Aqcd, and the collinear ones from nib down to y^mv^QCD- 
The distribution amplitude absorbs the soft contributions below the fac- 
torization scale (including the corresponding mass singularities). The 
hard-scattering kernel thus contains large logarithms coming from the region 
between m?, and a/?^V^^cd and it is these logarithms which we resum into 
the Wilson Coefficient C. The hard-scattering kernel has a factor of C, con- 
taining the resummed large logarithms and a remaining one-loop term Kr 
with no large logarithms. 
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fa) fb) 

Figure 10: Schematic representation of the large logarithms present in the 
contributions to the matrix element in the case of (a) the electromagnetic 
vertex and (b) the weak vertex. M represents the large scale, M ~ and 
Mb. 



The contribution from the weak current can then be factorized as follows: 



$ includes the soft logarithms of fijp/k"^. The hard-scattering kernel is the 
product of the coefficient C, which collects the resummed (exponentiated) 
logarithms of fip over large scales {rrib, q-), and of the "collinear" contribu- 
tions from the L^s, i.e. the logarithms of /i^^ ~ 2{q ■ k). If the factorization 
scale is set to = 0(mf,AQCD), the large logarithms of the matrix ele- 
ment are either included in the (non-perturbative) distribution amplitude or 
resummed in the coefficient C. 

Finally we combine the results from the corrections to the weak and elec- 
tromagnetic currents. We expect that the Wilson coefficient C, which resums 
the Sudakov logarithms, is associated with the weak current in any diagram 
and thus multiplies all the remaining contributions (the validity of this ex- 
pectation needs to be investigated at higher orders). We thus obtain the 
one-loop factorized expression for the form factors, including the resumma- 



l + ^^^^[L:' + Lf%k^;f.,) 
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tion of leading and next-to-leading logarithms: 



(83) 



where 



^, fBQuMBj_^ 

2V2E^ k+ 



1 + 



An 



and 



Kt{k+,EY, ftp) = log' 



2q_k^ 



TT 



^4) 



^5) 



Kt is the sum of two contributions: that from the electromagnetic current 
and Kj., the contribution from the weak current which has not been included 
in the resummation. 

In summary, following the formalism in ref. 



we have resummed the 

leading and next-to-leading logarithms of the form a" log"''"^(m;,/^/2A;+) and 
a" log"'(mfe/ \/2k+). Next-to-next-to- leading logarithms, of the form a"x 
log'^~^(mfe/\/2^+) have not been resummed here (although we do present the 
one- loop term of the form a^x constant in eqs. (|83|) - (|85|) ). These terms 



would require a three-loop calculation in the effective theory [jlT . 



6 Phenomenological Study 

In this section we briefly investigate the phenomenological consequence of 
the analysis presented above for the decay — >• e~^i"j (the discussion of the 
decay ii^v^ would be the same). The current experimental bounds 

on the branching ratios are T{B^ e^v))/TB < 2.0 x 10"'^ and T{B^ — * 
/i+z/7)/rij < 5.2 X 10-5 (both at 90% confidence level) pO[. is the total 
decay width of the 5-meson. Previous theoretical estimates of the branching 
ratios are typically in the range (1 — 5) x 10~^ [0, pl|] . 



The differential decay rate for the process B — >■ in terms of the form- 
factors Fv and Fa is given by [0: 



dEgdE^ 



aG||K.pM| 1^^^, ^ 2^^^ -R-1) + 1][F^ + Fl]{E,) 

-2R{1 -R){l + R- 2S)Fv{E^)Fa{E^)] , (86) 
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where R = 1 — 2E^/Mb, S = 2Ee/MB and and Ee are the energies of the 
photon and electron respectively. R and 5* satisfy < i? < 1 and i? < 5* < 1. 
Integrating over the electron's energy, we obtain the differential decay rate 



dV aGl\Vub\''M% 



dE^ 487r2 



R{1 - Rf[F^ + Fl]{E,) . (87) 



Our analysis of the form factors Fy and Fa-, is valid only if the energy of 
the photon is large. We therefore introduce a lower cutoff E'^ and consider 
the integrated decay rate 

/•Ms/2 

T{Ei)= / dE^ ^««^ 



j^c dEj 

E^ is formally of 0{Mb/2). 

At leading order of perturbation theory, the two form factors are given 

by 

^ 2E^Xb 1 — R Xb 
where is the first inverse moment of the 5- meson's distribution amplitude, 

' (90) 



Ab in 



Little is known about this quantity. Below, we will use the estimate of ref. 
\b = 350 ± 150 MeV Q. 

The integrated rate at leading order is 



where Rc = 1 — 2E^/Mb- Dividing the integrated rate by the experimentally 
measured total width (F^), we obtain the numerical LO estimate of the 
branching ratio {B^°{E^) = T^^{E^)/Tb) 

^ V3-6xl0-3y VlSOMeVy \ \b J 

xRl(l-^-^]. (92) 



®The definition of As here is the same as in ref. g, in spite of the different choice for 
the normahzation of 6^ . 
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The central values for \Vub\ and fb were taken from refs. [^] and ||22[ respec- 
tively. We therefore expect a fully integrated branching ratio F of a few 10~^, 
unless there is a significant enhancement of the soft-photon region of phase 
space (which is beyond the reach of our framework). 

It may be useful to normalize T{E^) by the purely leptonic decay rate 

r(i?+^/i+z/), 



/i+Z.) 



'PbMbM, 



in 



Ml 
Ml 



(93) 



since it absorbs the dependence on two relatively poorly known quantities: 
the S-meson decay constant, /b, and the modulus of the CKM matrix ele- 
ment \Vub\- (For the central values of \Vub\ and in eq. (p^) the branching 
fraction for the mode /xz/ is about 3.7 x 10~^. The current experimental upper 
bound is 2.1 X 10'^ at 90% confidence level |2|].) 

The integrated rate for the decay — > e^v^ is thus (at leading order 
and neglecting terms of 0{M'^/M1)) 



^(5^ 



aQl Ml Ml 
12vr Ml XI ^ 



0.22 



B 



A 



B 



1 - -i?, 

3 

1 - -R, 
3 



Ri 



3 



(94) 



(95) 



For the range of values of As given after eq. (|90|) , the product of the first two 
factors in eq. (^) varies in the range (25, 150). 

We now study how this result is modified at next-to-leading order in a^. 
It is convenient to rewrite the result for the form factors in eqs. (^) - (p5| ) 
in the form. 



Fa: 



IbQuMb 



47r 7 Xb 



Air 



X 



(1) 

B 



477 



A 



(2) 



(96) 



where a = — 7r^/6 — 1, 



log 



M, 



log(l - R) 



B 



(97) 
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and 

V2 rdk+^^,~ ^^_^V2KMb 



$^(Mlog" ^ : ■ (98) 



Xf Jo h ' f^l 

The integrated decay rate, up to next-to-leading order, is: 



2 



47r 47r At: 

where the A^^^'s are the ratios of inverse moments, A^^ = A^/A^''. C{jj,F) and 
L depend imphcitly on the integration variable R. 

We would like to estimate the size of the one-loop correction and the 
effect of re-summing the leading and next-to-leading logarithms. To this end 
it is necessary to re-expand the coefficient C{^f) up to one- loop order: 



^/ \ 1 OisCp 



2 log^ ^ - 4 log(l - i?) log ^ + 5 log 



Mb Mb Mb 

. (100) 



+21og2(l -R)+ 2Li2(i?) + ^^log(l -R) + ^ + Q 



R '12 

For any choice of input values for the cut-off Rc, the factorization scale fip, the 
large scale (which we have taken to be Mb here) and the non-perturbative 
inverse moments A^^ we can perform the integration in eq. (|99D numeri- 
cally. For illustration, in the following we choose = 1.25 GeV, and 
\b = 0.35 GeV (a change of the choice of the value for A^, for fixed A^i 
and N2, amounts to a rescaling of the decay rate in the discussion below). 

We now investigate how well we can evaluate the higher order corrections 
and how large they are. As a first estimate of the difference between the 
integrated decay rates calculated at LO and at NLO, we neglect the "small 
logarithms", x (k = 0,1,2). We will include these terms below, but 
it should be noted that in the summation formalism which we have used in 
sec. ^ they are of the same order as terms which have been neglected. From 
eq. (^), we see that now the NLO decay rate becomes independent of the 
inverse moments A^'' . The result is plotted in figure |Tl]. In this figure we have 
also included the one-loop result without resummation, replacing the Wilson 
coefficient C by its one-loop expression ( |100| ) in eq. ( P5| ) and neglecting all 
the terms of order x [and taking = OsIMb) in the resulting one-loop 
expression]. The figure shows the dependence of the integrated decay rate 
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20 



15 - 




Rc 



Figure 11: Integrated decay rate, normalized by the purely leptonic decay 
rate B fiu, as a function of the cut-off Rc = 1 — 2i?^, when we neglect the 
NLO contributions with small logarithms {as x L'^). The factorization scale 
is set to fip = 1-25 GeV and the large scale is taken to be Mb = 5.28 GeV. 
'LO" and 'NLO" denote the leading and next-to-leading results for the form 
factors, while 'NLO-fresum" denotes the result obtained after resumming the 
leading and next-to-leading logarithms. 



on the cut-off and demonstrates that the difference between the LO and (the 
resummed) NLO results is typically of the order of about 25%. 

For the remainder of this section we include the one-loop terms in the 
square parentheses in eq. (|99|) , and investigate the dependence of the decay 

(k) 

rate on the small logarithm L and on the inverse moments . Formally all 
the A^^ are of the same order as A^, and thus the A^^^'s are of 0(1). They can 
however, take significantly different values depending on the model used for 
the distribution amplitude. This is illustrated by the following two models 
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for the i?-meson distribution amplitude at the scale ij,f- 



1 

1 

B 
1 



A 



>(2) 



2k^ 

^0 

1 

a; 

1 

a; 
1 

a; 



■ exp 



^7 



log 



AqM 



B 



1 



log — 2 7 



TT 



(101) 

(102) 
(103) 

(104) 



and 



V2Kk: 



1 

A^ 
1 



A 



(1) 

B 
1 



(fc+ + i^'Ao)2 

1 

a; 

1 , V2XoMbK 

7- log 2 



A 



(2) 
B 



1 

a; 



For instance, taking As = Aq 
)wing values of the re 

I : Ni = -0.04, N2 = 1.81 



, 2V2X0MBK , 
log 2 + ^ 



0.35 GeV and K 



(105) 
(106) 
(107) 

(108) 
0.3 we obtain the 



following values of the ratios of inverse moments A^^ = \b/X^b'- 



II : Ni = -0.69 , N2 = 3.77 . (109) 



The first model was proposed in ref. 0, based on QCD sum rules. We 
propose the second one as an example of a distribution amplitude with well 
defined inverse moments, but divergent positive moments - as expected from 
the one-loop study of the renormalization properties of $^ above and in 
ref. m. 



In figure 12 we plot the difference (in percent) between the LO results 
and the NLO ones, as a function of the two ratios of inverse moments A^^i 
and N2 for —4 < < 4. Each band in the figure corresponds to a range 
of 20% and we present the values at the corners (iVi^2 = ± 4) in table |I|. The 
one-loop corrections are typically some (low) number of tens of percent and 
are seen to be sensitive to the values of the ratios of inverse moments, A^i,2- 
Such a significant variation of the rates with A^i_2 was to be expected from 
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(iVi,iV2) 


Rc = 0.4 
Resummed 1-loop only 


i?c = 0.6 
Resummed 1-loop only 


(-4,-4) 
(4,-4) 
(-4,4) 
(4,4) 


-24.6 % -16.8 % 
-52.7 % -35.0 % 
36.0 % 18.3 % 
-2.8 % -3.7 % 


-13.2 % -10.1 % 
-58.6 % -39.5 % 
51.0 % 26.2 % 
-11.7 % -9.1% 



Table 1: Percentage difference between LO and NLO integrated decay rates 
for various values of the two ratios of inverse moments A^^i and N2, for = 
0.4 and Rc = 0.6. In each case, we give the percentage of difference from 
the LO results in the case of the NLO resummed expression (|99|) , and the 
one-loop expression [with C expanded according to eq. (|100|) ]. 



an inspection of the integrand in eq. (^). The relative size of the one- loop 
corrections is also sensitive to the choice of the cut-off Rc- 

Even from our simple exploratory study it is clear that \Vub\fB/^B could 
only be extracted from the integrated rate of the radiative decay B '-/iu 
if we had precise information of the 5-meson's distribution amplitude (and 
in particular its shape) from an independent process. Such knowledge is 
necessary to estimate the potentially significant contribution of the inverse 
moments A^'* and A^'' to the NLO prediction. 

We should stress again that the resummed formulae in eqs. (p6D and (|99|) 
are not complete in that we have not resummed the NNLO logarithms, which 
at one-loop order are the terms which are proportional to ag without any large 
logarithms. 

An interesting question, but one which is beyond the scope of the present 
study, is whether it would be possible to invert the above procedure and 
determine the important features of the 5-meson's distribution amplitude 
from the measured (in future) decay distribution with sufficient precision for 
use in calculations of decay rates of other processes. The B — > jiui radiative 
decay, with no hadrons other than the 5-meson, seems to be a particularly 
appropriate choice for such a determination of the distribution amplitude. 

7 Conclusions 

In this paper, we have studied the radiative decay B jiui at one-loop 
order in the framework of QCD factorization. We have explicitly verified 
that factorization is valid at this order and in eqs. ( [55| ) and (^) we present 
the result for the hard-scattering amplitude. The non-perturbative physics is 
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contained in the 5-meson's distribution amplitude, which is defined on the 
hght-cone (i.e. z+ and zx_ are set to zero in eq. (0) ) and, at least at this order 
of perturbation theory, there is no need to re-formulate the factorization 
formalism in terms of wave functions which depend also on the transverse 
momenta. On this point we disagree with earlier claims |0] . 

In ref. the perturbative expansion of the light-cone Bethe Salpeter wave 
function is not included in the matching procedure and in addition, the wave 
function is redefined by a perturbative factor k±) in order to modify 

the evolution behaviour of the wave function (see eq.(26) in ref. 0); the hard 
scattering kernel therefore is redefined by a factor of 1//. This corresponds to 
trying to absorb (at least part of) the contributions from $(^) '^°^(g)T'^*') (which 
depends on the external transverse momentum) and <|)(^)™'^(g)T*^°) into T(i) by 
redefining the wave function. However, we stress that such a redefinition is 
performed perturbatively and that the non-perturbative physics is contained 
in the distribution amplitude (0). 

The renormalization-group properties of the distribution amplitude of 
the 5-meson are very different from those of light mesons and cannot be 
obtained from the positive moments Q. This makes the determination of 
the distribution amplitude using standard non-perturbative techniques, such 
as lattice simulations or QCD sum rules, considerably more difficult. On 
the other-hand, one should explore further the extent to which one might 
use future experimental data on the photon energy distribution m B ^ 'jii'e 
radiative decays to determine the properties of the distribution amplitude 
which are needed in the phenomenology of B decays to two light-mesons. 

The resulting one-loop hard-scattering kernel contains large double and 
single logarithms of the ratio mh/k+, due to Sudakov effects at the weak 
b —* u vertex. These logarithms are precisely those which one obtains from 
the Soft-Collinear Effective Theory (SCET) ^ (see refs. 0, 0, |2| and 



references therein for discussions of the resummation of Sudakov logarithms 
using the Wilson Line formalism). Assuming that this is also the case at 
higher orders, we can resum the large logarithms using the formalism devel- 
oped in ref. [O]. Moreover, we have been able to match the SCET diagrams 



involving soft gluons with the contributions to the distribution amplitude. 
Our phenomenological study shows that the corrections are expected to be 
significant, up to 50% or so depending on the values of the parameters A^-* 
defined in eq. (|98|) . Indeed, as mentioned above, one can envisage in princi- 
ple using B — i> 'jii'i decays to determine these parameters and to use them 
in calculating predictions for other processes, such as two-body non-leptonic 
-B-decays. 

An important motivation for this study is the need to develop the QCD 
factorization formalism for two-body non-leptonic 5-decays. In particular. 
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it is necessary to understand the hard spectator interactions (the term on 
the second hne of eq. (|lD ) beyond the tree-level, and the corresponding loop 
corrections share the key features of the calculations described in this paper. 
These include the dependence on the three scales and the need to control the 
Sudakov large logarithms. 

Our conclusions for the development and exploitation of the QCD Factor- 
ization formalism are very optimistic but a considerable amount of work still 
needs to be done. The explicit demonstration of factorization at one-loop 
order in this paper, needs to be extended to higher orders. Most probably 
techniques such as those incorporated into the SCET will be very useful in 
this context. These studies have also to be applied to the decays of 5-mesons 
into two light mesons, so that the large (and growing) amount of experimen- 
tal data, predominantly from the 5-factories, can be analyzed in terms of the 
fundamental parameters of QCD (in particular the CKM-matrix elements) 
and provide an understanding of CP-violation in the quark sector. 
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2 4 



Figure 12: Difference between LO and NLO integrated decay rates as a 
function of the two ratios of inverse moments Ni and N2, for Rc = 0.4 
(upper line) and Rc = 0.6 (lower line). In each case, the plot on the left 
hand-side corresponds to the NLO resummed expression (PPj), and the plot 
on the right hand-side to the one-loop expression [with C expanded according 
to eq. QTUDI)]. 
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